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Abstract 

We study the tunneling of Dirac fermions in graphene through a double barrier potential allowing 
the carriers to have an effective mass inside the barrier as generated by a lattice miss-match with the 
boron nitride substrate. The consequences of this gap opening on the transmission are investigated. 
The realization of resonant tunneling conditions is also analyzed. 
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1 Introduction 



After the experimental realization of graphene in 2005, this system became an attractive subject not 
only to experimentalists but also to theoretical physicists. This is partially due to the relativistic-like 
behavior of its massless carriers. In addition to its anomalous quantum Hall effect [UE]; graphene pro- 
vides an example of condensed matter systems where quantum electrodynamics tools can be applied [3] . 
These new developments offered a condensed matter laboratory to perform many investigations, which 
were solely available to high energy physicist in the past. Because of this similarity with relativistic 
fermions, graphene systems triggered the interest of the scientific community for studying massless 
Dirac fermions in two-dimensional spaces. 

One of the characteristics of Dirac fermions in graphene is their ability to tunnel through very high 
potential barriers with unit probability [H [5], contrary to the usual intuition. This so-called Klein 
tunneling of chiral particles has long ago been proposed in the framework of quantum electrodynam- 
ics O El [H], but was never observed experimentally. However, as appealing as the Klein tunneling 
may sound from the fundamental research point of view, its presence in graphene is unwanted when 
it comes to applications of graphene because space confinements of the carriers is of great importance 
in nanoelectronics. One way to overcome these difficulties is through generating a gap in the energy 
spectrum or equivalently to generate a mass term in the associated Dirac equation. 

Experimentally it has been shown that one can generate a gap in the graphene spectrum if a 
graphene sheet is deposited on top of hexagonal boron nitride (BN) [9]. In graphene the carbon- 
carbon distance is 1.42 A, when deposited on a substrate of boron nitride, a band gap insulator with 
a boron to nitrogen distance of the order of 1.45 A [lOj creates a lattice miss- match which is at the 
origin of gap opening in the graphene sheet. It was shown that the most stable configuration obtains 
when one carbon is on top of a boron and the other carbon in the unit cell is centered above a BN 
ring, the value of the induced gap is of the order of 53 meV. Depositing graphene on a metal surface 
with a BN buffer layer leads to n-doped graphene with an energy gap of 0.5 eV |11] . 

Theoretically, tunneling phenomena through graphene deposited on Si02-BN substrate in zero 
magnetic field was discussed |12j . In that work it was assumed that it was possible to manufacture a 
thin layer with Si02-BN interfaces, on top of which a graphene flake was deposited. This arrangement 
induces spatial regions where graphene has a vanishing gap intercalated with regions where BN will 
cause a finite gap. The graphene physics was considered in two different regions: the /c— region, where 
the graphene sheet is standing on top of Si02, and a g— region where the sheet stands over a BN sheet 
hence causing a mass-like term and inducing an energy gap of value 2t' . The effect of chiral electrons 
in graphene through these region was studied, the electronic spectrum changes from the usual linear 
dispersion to a hyperbolic dispersion, due to the presence of a gap. It was shown that contrary to the 
tunneling through a potential barrier, the transmission of electrons was, in this case, smaller than the 
one for normal incidence. These results were extended to the case of a magnetic field case |21J and 
different generalizations were obtained. 

Motivated by the above reasons and in particular the investigation made in |12j and |21j . we would 
like to extend the study of tunneling phenomena through a graphene sheet deposited on Si02-BN layer 
in the presence of an external magnetic field B and double barrier potentials. Note that, actually there 
were many theoretical investigations on tunneling of electrons in graphene in the presence of magnetic 
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field e.g. [131 IHl US ES]. However, they did not consider a graphene system deposited on Si02-BN 
and hence their graphene system was gap less. 

More precisely, we consider a system made of two different regions, where the second is charac- 
terized by an energy gap t', in a perpendicular magnetic field. The energy spectra are obtained in 
both regions in terms of two different Landau levels and t' . For concreteness, we consider a barrier 
in magnetic field and study the tunneling effect. Indeed, from the continuity condition we obtain dif- 
ferent solutions, which allowed us to explicitly determine the reflection and transmission coefficients 
in different regions those are then shown to obey the probability conservation law. We study the 
transmission behavior in terms of the energy ratio ^ and discuss the possibility of resonant tunneling. 

The present paper is organized as follows. In section 2, we formulate our problem to include 
different part in the Hamiltonian describing the system made of graphene. The corresponding solution 
of the energy spectrum will be the subject of section 3. The scattering problem for Dirac fermions will 
be solved in section 4 where the transfer matrix will be explicitly given. In section 5, we determine the 
reflection and transmission coefficients using the associated current density. We conclude our work in 
the flnal section. 



2 Problem formulation 

To deal with our task let us flrst describe the geometry of our system, which is made of flve regions 
denoted by j = 1, 2, • • • 5. Each region is characterized by its potential and interaction with external 
sources. The potential proflle is summarized in Figure 1: 



V 



U 
B 

t' 



-do 



-d\ di 



dj_ 



Figure 1: The energy diagram of the potential V{x) describing our double barrier along with its physical 
parameters. 



The double barrier potential V{x) is defined by 



V{x) 



0, 



di < \x\ < d2 
\x\ < di 
otherwise. 



According the region labels, we denote the corresponding constant potential values by Vj. We introduce 
a uniform perpendicular magnetic field, along the z-direction, constrained to the well region between 
the two barriers and is defined by 



S(x,y) = SoG [dl-x^) 
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where Bq is the strength of the magnetic field within the strip located in the region |x| < di and B = 
otherwise, G is the Heaviside step function. Choosing the Landau gauge and imposing continuity of 
the vector potential at the boundary to avoid nonphysical effects, we end up with the following vector 
potential 

— di, X < —di 

X, \x\ < di (3) 



_ di, X > di 



with the magnetic length defined hy Ib = \/ c/eBo, c being the speed of light and e the electronic 
charge. 

Once we have defined the potential parameters relevant to all regions, we introduce the Hamiltonian 
for one-pseudospin component describing our system. In region j = 1, 2, 4, 5 it can be written as 

Hj = Vfa -TT + Vj (4) 

while in region (3) we have 

H3 = vfza -TT + U + ta^ (5) 

where is the Fermi velocity, tt = p + eA/c is the two-component momentum, t' is the mass term, 
<? = {o'x,cry) and fi^ are the usual Pauli matrices. To proceed further, we need to find the solutions of 
the corresponding Dirac equation and therefore the energy spectrum. 

3 Energy spectrum 

We determine the eigenvalues and eigenspinors corresponding to the Hamiltonian Hj in different 
regions. Indeed, from dH) we can write in regions (1, 2, 4, 5) the Hamiltonian as 



Vj/vf pjx -ipy + i^ 



where the parameter d is defined by 



d=H^' (7) 
I — di, X > di. 

To solve the eigenvalue problem, we can separate variables and write the eigenspinors as 

i^j{x,y) = ipj{x)x{y) (8) 

where ^Pj{x) = fj-Y: 'Pj± are the upper and lower spinor components, and xiu) = e*'^^^, with ky 

a real parameter that stands for the wave number of the excitations along the y-axis. After rescaling 
energy e = E/vf and potential Vj = Vj/vf, the resulting reduced time independent Dirac equation is 
given by 



Pjx - iPy + ijT 




Pjx + iPy-ijT j \ ipj^ J \ ifj. 



(9) 
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This eigenproblem can be written as two linear differential equations of the from 

Pjx - iPy + i-^^ fj- = (e - Vj)(pj+ (10) 

(pjx + iPy - = (e - Vj)<fj- (11) 

where pjx is the x-component of the momentum operator in the j-th potential region. The above 
equation can be solved, using the unit system (h = c = e = 1), to get the energy eigenvalues 



e-vj = sj\{kjxf+[ky--^] (12) 



where Sj = sign(e — Vj) is the usual sign function, equal to ±1 for a positive and negative argument, 
respectively. It is clear that the wave vector along the x-direction can be written as 



= y (e - VjY - (^ky - . (13) 

Consider now a particle-like scattering state (e > 0), the particle is incident from the left side in 
the j-th potential region, with incoming wave vector kj = {kjx,ky) and position r = (x,y). Thus, the 
incoming spinor in the j-th potential region is given by 



I'll Jk,.r 



where Zj is given by 



]J {k,.y + {ky - 



Sjc'^^ (15) 



d 



and the phase is (pj = arctan I ^ ^ I . The sign sj = ±1 again refers to conduction and valence 



bands in region j = 1,2,4,5. The corresponding left and right moving spinors are defined by 

Jikj^x+kyy) ^26) 















^1 





iPj-{^,y) = _ j e^^'^.^^+M). (17) 

To complete the solutions of the energy spectrum, we consider now region 3 where |x| < di, 
T^x = Px, T^y = Py + -^yi ^^id vector potential is A{x) = x//^, then ^ becomes 

/ U/v^ + t/vf -idx-iky-i^\ 

H^=\ , ^ . (18) 

\-idx + iky + i^ U/v^-t/v^ j 

To diagonalize this Hamiltonian we proceed by defining the usual bosonic operators 

^(dx + ky + ^], a} = ^( ' ^ 



a = ^[dx + ky + -j^y af = — y-dx + ky + -j^] (19) 
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which satisfy the commutation relation [a, a^] = 1. RcscaHng our energies fj, = t /vp and u = U/vf, 
we can write the Hamiltonian in terms of a and as follows 

(u + a —i^a \ , , 

The eigenvalue equation for the spinor ip3{x,y), that is H^tp^lxjy) = Eip3{x,y) gives Jp3(x,y) = 
f3ix)x{y) with x{y) = a"-^"^ and Lpz{x) = (^33+ v?3_)* verifies 



's c 



(21) 



where again e = E/vf. The eigenproblem can then be rewritten as two relations between the two 
spinor components 

V2 

{u + n)ip3+ - i—aip3- = eip3+ (22) 

IB 

i—a^ip3+ + {u- I^)ip3- = eip3- (23) 

IB 

which lead to an eigenvalue equation for the operator aa} and (^3+ 

2 

[(e - uf - /x^] (p3+ = -^aa^ip3+. (24) 

It is clear that (p3+ is an eigenstate of the number operator N = a^a and therefore we identify (^53+ 
with the eigenstates of the harmonic oscillator |n — 1), namely 

ip3+^\n-l) (25) 

and the corresponding energy spectrum reads 

{e-u) = ±en = ±^^/{id^)^+2^. (26) 
The second spinor component can now be written as follows 

iV2a^ I i\ I \ /o'7\ 

V>3- = 7 n r~ ~ 1/ = 7 r; r~ ) 

[e — u)Ib + IJ-Ib [e — u)Ib + IJ-Ib 



where ■\/2n = ^/{enlBY ~ (m^b)^! which gives 



<^3- = W^I")- (28) 

Where the it signs correspond to positive and negative energy solutions, respectively. After normal- 
ization we obtain the eigenspinors for positive and negative energies 



1 / , 
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By introducing the parabolic cylinder functions Dn{x) 
barrier region | a; |< di can be expressed as 



the solution in the 



1 

72 



± V2 + kyls 



Jkyy 



(30) 



We should note that the above it signs stand for waves traveling to right and left, respectively, along 
our usual convention for traveling waves g^*'^-''^^. In this context the full solution in region 3 is given 
by a linear combination, c+ times ()30p with the upper sign added to c_ times ()30p with the lower sign 
where c± are arbitrary coefficients to be fixed by the boundqry conditions as exposed in the following 
section. The above results (|16p -(|17 p and (j30p summarizes the solutions of the energy spectrum in 
different regions composing our graphene sheet. Next, we will show how to implement the results 
obtained so far to study the tunneling of Dirac fermions through a double barrier potential. 



4 Boundary conditions 

It is straightforward to solve the scattering problem for Dirac fermions. We assume that the incident 
wave propagates at an angle (pi with respect to the x-axis, where the superscript indicate the potential 
region. Let us first write explicitly the eigenvalues and the corresponding eigenspinors with reflected 
and transmitted parts in each region. In the incident region for x < —d2 (rigion 1), we have 



\h,r+{ky-^X (31) 



V2\ziJ V2\-zlJ 

where r is the reflection amplitude. It is clear that the shift in ky is due to our gauge choice for the 
vector potential. It is convenient to parameterize the momenta as follows 

fell; = e cos (/>i , ky = e sin (pi + . (33) 

'■B 

In the barrier region —d2 < x < —di (region 2), we have 

'2 



e-v = s^^j {k2xY + [ky 



di 



(34) 



^2 = — ( ''^ I 6**^*^2^^"'''^^^^ + — ( ^ I e^^~^'^'^^^^'"'^^ (35) 
V2\ Z2 J V2\-z^J 

where a and b are two wavefunction parameters and v = V jvp. The corresponding, suitably parame- 
terized, momentum is given by 

fc2x = (e - t^) cos(/)2, fej/ = (e - v) sin(/)2 + (36) 

^B 



6 



Similarly for the barrier region di < x < d2 (region 4) we have 



e-v = S4\l {k2xY + ( ^2/ + ^ 



-04 



V2 



1 



,i{k4xX+kyy) 



+ 



V2 I -zl 



J{-k4xX+kyy) 



(37) 
(38) 



where e and / are the corresponding wavefunction parameters. Similarly, the momentum can be 
defined by 

k^x = (c ~ ^) cos (j)4, ky = {e — v) sin ^4 — (39) 

Finally, in the transmission region, x > d2 (region 5), we have only transmitted waves 

1. 



e= \j {k^xf + {ky + 



1 

^5 



J{kzxX+kyy) 



where t is the transmission amplitude, the corresponding momentum reads as 

k^x = e cos <^5, ky = e sin <^5 - ^ • 



(40) 
(41) 

(42) 



The refraction angles 02, ^4 and ^5 at the interfaces are obtained by requiring conservation of the 
momentum py. This leads to a simplified expression of these angles in terms of (f)\ 



sm 02 
sin 04 



e — V 
e 

e — V 



■ sm 01 



■ sin 01 + 



2di 



(e - v)l 



•A • ^ ^ 2di 

sm 05 = sm 01 + . 

€1 r. 



(43) 
(44) 

(45) 



We should point out at this stage that we were unfortunately forced to adopt a somehow cumbersome 
notation for our wavefunction parameters in different potential regions due to the relatively large 
number of necessary subscripts and superscripts. 

Before matching the eigenspinors at the boundaries, let us define the following shorthand notation 
for the dimensionless parameters 



0'n± 



V2 



1±-, dn±-- , 



and the parabolic cylindric functions by 

??i± = D^{euiBr-{i^iBn/2-i 



±V2 



-di 



6± = -D[(,^,^)2_(^,^)2]/2 



V Is 
di 



ylB 



±V2 {-j^ + kyh 



(46) 

(47) 
(48) 
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Now, requiring the continuity of the spinor wavefunctions at each junction interface give rise to a set 
of equations. Indeed, at the point x = —d2- 



X = —di: 



X = di: 



X = do: 



^—ikixd2 _j_ j,^ikixd2 _ g^^—ik2xd2 _j_ Jjgik2xd2 
^^-ik2xdi _|_ 5g«A:2^di ^ c+an+r]i + C-an-Tj'^ 

{ggik4xd2 _j_ Jg—ik4xd2 _ ^gjA;5a;l^2 
g^^gjfc4.rf2 _ f^*^-ik4xd2 = ^^j.gifc5.rf2 



(49) 



(50) 



(51) 



(52) 

In the above formulae, the parameters ri2± and ^2± si's defined similarly to ()47p -()48p but with di 
replaced by —di. All these equations can be written in compact form by introducing the transfer 
matrix M, which is given by 

/ _^ik\xd2 Q-ik2xd2 ^ik2xd2 \ 

±f,iklxd2 I ZLLeik2xd2 

Zl Z2 

e-ifc2.di gifc2.di -a„r?i+ -a„_?]i_ 

^2e~*''2.di =±eik2xdi -id^^^^^ idnM- 

-an+m+ -an-V2~ e'^*^'^^ (.-ikixdi 

idn+i2+ -idn-^2- -z^e^P^-'^^ _Le-ifc4.di g 

Q Q Q gik4,xd2 Q—ik4,xd2 _^ik5xd2 

y Z^e'''^^'^^ _Xf.-ik4xd2 _^^^iksxd2 J 

(53) 

Defining $ = (r, a, 5, c+, c_, e, /, t)* as a vector quantity made of all unknown wavefunction param- 
eters in our problem and H = (e"*'^^^'^^, zie"*'^^"''^^ , 0, 0, 0, 0, 0, 0)*, where the superscript t stands for 
transpose, the compact form of equation ()49p - (j52p can then be written as 

M^> = H. (54) 

Of course the same problem could be formulated in terms of 2 x 2 matrices if we use the concept 
of transfer matrix from one potential region to another. This last formulation will be much more 
adequate in dealing with periodic systems and applying the Bloch theorem to find the associated 
energy bands. 



5 Reflection and transmission coefficients 

Now we are ready for the computation of the reflection R and transmission T coefficients. For this 
purpose, we introduce the associated current density J, which defines R and T as 

i? = ^, T=^-^ (55) 
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where Jin, Jre and Jtr stand for the incident, reflected and transmitted components of the current 
density, respectively. It is easy to show that the current density J is given by 

J = ev^ij)'^ Gxi^ (56) 
which gives the fohowing results for the incident, reflected and transmitted components 

Jm = ev^{lplin)^ CTa^lpiin = CVf . =^ (57) 

^(M^+(fc,-|) 

Jre = ev^{il)ire)^axipire = ev^r*r—= ^^"^ ^ (58) 

^(M^+(fc,-l) 

Jtr = ev^{ijj^tT)^ax'4'ur = evft*t—= ^^"^ (59) 

^ (fc5,)2 +[ky + j 

Injecting these results in (|55p we obtain 



ik,x? + - ^) 



2 



t|2 (60) 



klx LZ, /, d,N2 



R = |r|2. (61) 



Now using the conservation of energy 



"1 \ hi \o . I , "1 ^ 



{klxY +[ky-fj = y {k,x? +[ky + ^^] (62) 



we find the constraint 



klx = J {k5x)^ + 4ky^ (63) 



"B 

which allows us to express the transmission probability T in the following form 



k5x |ji2 



A ■ (64) 



'{k^xr + Aky^ 

To get an explicit expression for T, we should determine the transmission amplitude t. After some 
lengthy algebra, one can solve the linear system of equations ()54p to obtain the transmission and 
reflection amplitudes in closed form 

iV2 e^('^l+'i2)(*^2.+fc4.)(l + (zi)2)(i + (Z2)2)(l + (Z4)2) 



Isle - u\ e^d2(k^_^+k5,:) [^(1 + Z1Z2) + B{zi - Z2)] 

g-2id2fci.^^ [^(^1 - Z2) - B{1 + Z1Z2)] 



(%+6- + m-i2+) (65) 

(66) 



^(1 + zi^s) + -8(2:1 - Z2) 
where we have defined the following quantities 

A = e2^('^i^2.+rf2fc4.)(^^_^^)(„ + ^^2^4 + 724-^2) (67) 

+g2*di(fc2.+P4.)(i + ^^^^) + /3^2 + 7 + '522 24) 

S = e2^'^2(fc2.+fc4.)(^^_2^)(_„^2 + /3z4-7Z2Z4-<5) (68) 

+g2i{d2fc2,+difc4.)(i + (^^2^4 + /3 - 7Z2 + fc4) 
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as well as 



e — n + /i 

a = (r/i_?72+ - ??i+??2-) (69) 

e — u 

P = r-( w —^(6-6+ -6+6-) (70) 

7 = r(m-6+ + m+6-) (7i) 

Ib\€ — U\ 

■ (7/2-6+ + r/2+6-)- (72) 



/file — u 



The resulting reflection and transmission coefficients have been computed numerically and are plotted 
in Figures 3, 4, 5, 6 for several parameter values (e, u, di, (i2, [x)- For instance a typical value of 
the magnetic field say = 4T, the magnetic length is Ib = 13 nm, and els = 1 corresponds to the 
energy E = 4A meV [I3]. 

Now, let us study the reflection and transmission coefficients versus the energy elB- The quantity 
kylB = rn* plays a very important role in the transmission of Dirac fermions via the obstacles created 
by the series of scattering potentials, because it associates an effective mass to the particle and hence 
determines the threshold for the allowed energies. However, the application of the magnetic field in 
the intermediate zone seems to reduce this effective mass to (^kylB — in the incidence region while 

it increases it to (^kylB + in the transmission region as shown in Figure 2. The allowed energies 
are then determined by the greater effective mass, namely gIb ^ kylB + j^- 















h + — ^ 






^ 4 





Figure 2: The energy configuration of the double barrier potential illustrating the effect of inhomogeneous 
effective mass in the incidence and reflection regions. 



If we cancel the well region by setting di = then we reproduce the usual single barrier transmis- 
sion as reflected in figure 3 which shows clearly the Klein zone followed by a full reflection zone. On 
the other hand if we keep the well region and cancel both the applied magnetic field and mass term 
in the well region, the series of potentials behaves like a simple double barrier with the same effective 
mass ky all over. Thus, in this case we reproduce exactly the transmission obtained in [22], for the 
massive ID Dirac equation with m = ky. 
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T 




Figure 3: Transmission and reflection coefficient as a function of energy cIb for a single barrier with ^ = 0, 
^2^5 = 1-5, vIb = 6 and kyls = 1 

The result shown in Figure 3 can be interpreted as follows. The transmission as function of energy, 
through a single barrier, starts from the minimum required energy equal to the effective mass m* = 
kylB then oscillates reaching full transmission due to constructive interference in the (Klein zone) for 
energies between m* and V — m* . This is then followed by a gap region between V — m* and V + m* 
where the transmission drops to zero. For energies above V + m* the transmission has the usual 
oscillating behavior then reaches full total transmission asymptotically. 

In Figure 4 we show the influence of the size of the well region and magnetic field on the trans- 
mission. First we consider the same conditions as above but we vary di to obtain Figure 4a: 



T 




Figure 4a: Transmission as a function of energy elB with = {0.05, 0.1, 0.2}, ^ = 1.5, VIb = 6, 

uIb = 4, filB = 4, kylB = 1. 

According to Figure 4, the size of the intermediate zone affects the form of transmission in these four 
zones. As di increases, the effective mass increases and the oscillations in the Klein zone get reduced. 
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This strong reduction in the transmission in the Klein zone seem to suggest the potential suppression 
of the Klein tunneling as we increase di . On the other hand the size of the gap region increases as we 
increase di. However, this is not the case as we increase ^Ib where there is some change amplitude as 
shown in Figure 4b: 




It is worthwhile to analyze the reflection and transmission coefficients versus the potential. In 
doing so, we fix the energy and choose a value of di, then we compute the transmission as shown in 
Figure 5a and 5b: 
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Figure 5b: Transmission as a function of potential uIb for = {0, 0.1, 0.2, 0.3, 0.5, 0.6, 0.7} with 
= 1.5, cIb = 6, vIb = 4, fils = 4 and kylB = 1. 

It is clearly seen that for a given energy, the transmission decreases if di increases and then vanishes. 

Let us examine the transmission coefficient as function of angle incidence This is shown in 
Figure 6: 




Figure 6a, 6b, 6c: Polar plot of a curve with radius (transmission T) as a function of angle (pi with: 

(a) : g- = {0.5, 1.5, 3, 3.67} with eh = 3.7, = f^' = uIb = 0, fils = 0, kylB = 1. 

(b) : elB = {1, 2, 4, 8} with = 0.5, g = 0.6, vIb = 0.5, uIb = 0.4, ^Ib = 1, kyh = 1. 

(c) : g = {0.5, 1.5, 3, 3.67} with eh = 2, g = 0.6, vIb = 0.5, uIb = 0.4, fih = 1, Kyh = 1. 

Figure 6a reproduces exactly the result of De Martino et al. [13j, this reference was the first to treat 
the confinement of Dirac fermions by an inhomogeneous magnetic field. This polar graph shows the 
transmission as a function of the incidence angle, the outermost circle corresponds to full transmission, 
T = 1, while the origin of this plot represents zero transmission. We note that for certain incidence 
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angles the transmission is not allowed, in fact for els < kyls + all waves are completely reflected. 
It is worth mentioning that the transmission is uniquely defined by the incidence angle i.e. each radial 
line representing a given incidence angle intersects the transmission curve at one point. In Figure 6a 
we see that the transmission vanishes for values of below the critical value els = 6.7. In Figure 6b 
we fix ^ and ^ and vary the energy, we observe that the transmission get reduced as we increase the 
increase the energy. In Figure 6c we fix ^ and vary f^, the transmission vanishes below the critical 
value = 3.7 for our parameters. 




Figure 6d,6e: Polar plot of a curve with radius (transmission T) as a function of angle 0i with : 

(d) : g = {0.6, 0.1, 2, 3} with eh = 8. = 0.5, vIb = 0.5, uh = 0.4, fih = 1, Kyh = 1. 

(e) : filB = {0.5, 1, 2, 3} with els = 4, ^ = 0.5, J = 0.6, vIb = 0.5, uIb = 0.4, Kyh = 1. 

In Figure 6d we fix and vary p-, the transmission is less sensitive to the variations of and shows 
an oscillatory behavior which becomes more apparent for larger values of j^. Finally in Figure 6e 
we show how the transmission is affected by the effective mass term reflected by /x, the transmission 
decreases sharply as we increase /i. 



6 Conclusion 

To conclude, we have studied the effect of gap opening on the transport properties of graphene Dirac 
fermions. The effective mass is generated by a lattice miss-match between the boron nitride substrate 
and the 2D graphene sheet. An additional homogeneous magnetic field localized within the well 
potential region causes the effective mass to change and affects asymmetrically the incidence and 
reflection regions. 

We have found that in contrast to electrostatic barriers, magnetic barriers are able confine Dirac 
fermions, a property of great importance from the practical point of view. In fact as can be seen 
in Figure 6, the transmission reduces a lot as we increase the effective mass parameter, t' in our 
original model (5), or the width, di, of the effective mass region. These results were exposed clearly 
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on the polar graph showing the transmission as a function of the incidence angle in Figure 6. The 
outermost semicircle corresponds to unit transmission is not intersected by the transmission curve, 
even at normal incidence, if we increase the effective mass or width of the magnetic field region. 

We hope that the present interplay between gap opening and local magnetic field effect will give 
more freedom to experimentalists to develop graphene based electronic devices. 
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